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1. INTRODUCTION 

Consider the focusing cubic nonlinear Schrodinger (NLS) equation 

iut + Uxx + 2\ufu = 0, (1.1) 

where subscripts denote the appropriate partial derivatives. It arises in applications as 
diverse as wave propagation in nonlinear media [15], surface waves on suflBciently deep 
waters [14,15], and signal propagation in optical fibers [10-12]. Its initial-value problem 
is known to be solvable by the inverse scattering transform method [1,2,13,15]. In other 
words, certain solutions to (1.1) can be viewed as a potential in the Zakharov-Shabat 
system 

<fi{X,x,t), (1.2) 

where an asterisk is used to denote complex conjugation, and u{x, t) can be recovered from 
u{x, 0) with the help of the scattering data sets for (1.2) at t = and at time t. 

Exact solutions to nonlinear partial differential equations are of great interest. Such 
solutions may be helpful to better understand the corresponding nonlinearity, and they may 
also be useful in producing testing means to determine accuracy of numerical methods for 
solving nonlinear partial differential equations. This paper is related to exact solutions to 
(1.1). 

In previous papers [3,6-8] we presented a method to construct exact solutions to (1.1) 
that are globally analytic on the entire xt-plane and decay exponentially as a; — > ±oo at 
each fixed t e R. This has been achieved by using a matrix triplet {A,B,C), where all 
eigenvalues of A have positive real parts. A similar method was applied to the Korteweg- 
de Vries equation on the half-line [4]. The same method is also applicable to various 
other nonlinear partial differential equations that are integrable by the inverse scattering 
transform with the help of a Marchenko integral equation. 

In this paper we analyze the method of [3] when the eigenvalues of the matrix A in 
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d(f{X, X, t) 
dx 



—iX u{x,t) 
—u{x,t)* iX 



the triplet {A,B,C) do not all have positive real parts. It is already know that if one or 
more eigenvalues of A are purely imaginary, the corresponding scattering coefficients for 
(1.2) contain discontinuities at some real values of A and hence the corresponding u{x,t) 
cannot be analytic on the entire a;t-plane. Furthermore, it is already known that for soliton 
solutions to (1.1), the corresponding transmission coefficients for (1.2) must have a pole 
and a zero appearing as a pair located symmetrically with respect to the real axis, which 
implies that a pair (or more pairs) of eigenvalues of A cannot be located symmetrically with 
respect to the imaginary axis. Thus, in our paper we mainly concentrate on the case where 
eigenvalues of A may occur anywhere on the complex plane, but no eigenvalues of A are on 
the imaginary axis and no pairs of eigenvalues are located symmetrically with respect to 
the imaginary axis. For such triplets {A, B, C) we show that there is an equivalent triplet 
{A, B, C) yielding the same solution u{x, t) to (1.1), where all eigenvalues of A have positive 
real parts. The corresponding solutions u{x,t) are then analytic on the entire a;t-plane 
and they are soliton solutions with any number of poles in the corresponding transmission 
coefficients and with any multiplicities of such poles. For such triplets we also explicitly 
evaluate the corresponding transmission coefficients, bound-state norming constants, and 
the corresponding Jost solutions to (1.2). We also consider the generalization of our results 
to the matrix case, where the scalar quantity u{x, t) in the NLS equation is replaced by a 
matrix-valued function of x and t. 

Our paper is organized as follows. In Section 2 we present the preliminary material 
by providing an outline of the method of [3] and introduce exact solutions u{x,t) to (1.1) 
constructed via a triplet {A,B,C). In Section 3 we exploit a certain symmetry in such 
exact solutions and show that some (or all) eigenvalues of A can be chosen either on the 
right or on the left half complex plane without changing u{x,t). In Section 4 we show that 
such solutions are solitons with any number of poles in the corresponding transmission 
coefficients and with any multiplicities of such poles. We also explicitly evaluate the 
corresponding transmission coefficients and bound-state norming constants and the Jost 
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solutions to (1.2). In Section 5 we show that the results of Sections 2-4 obtained for the 
(scalar) NLS equation (1.1) remain valid for the matrix NLS equation 

iut + Uxx + '^uu^u = 0, (1-3) 

where u is now m x n matrix valued, the dagger denotes the matrix adjoint (matrix 
transpose and complex conjugate), and the associated Zakharov-Shabat system is given 

by 

dif{X,X,t) \ -i^Im U(x,t)] 

ip{X,x,t), (1.4) 



—iXIjn u(x, t) 
-u{x,ty iXIn 



dx 

with /„ denoting the nxn identity matrix. Finally, we conclude in Section 6 with an explicit 
example, which was earlier studied as Example 7.2 of [3]. Since one of the eigenvalues of 
A in that example has a negative real part, we earlier conjectured that it might be a 
nonsoliton solution. Using our current results, we verify in Section 6 that it is a two- 
soliton solution and we explicitly evaluate the corresponding transmission coefficients, the 
bound-state norming constants, and the Jost solutions to (1.2). 

Note that (1.1) and (1.3) are identical when tt is a scalar, and similarly (1.2) and (1.4) 
are identical when m is a scalar. Throughout our paper we refer to (1.3) as the scalar NLS 
equation when u is a scalar and as the matrix NLS equation when u is matrix valued. This 
convenience enables us to state all our results for the NLS equation so that, as shown in 
Section 5, they remain valid whether the scalar case or the matrix case is considered. 

2. PRELIMINARIES 

In this section we establish our notation and provide the preliminaries for certain 
exact solutions to the focusing NLS equation. Such exact solutions are expressed [3] in 
terms of a triplet of matrices {A, B,C). 

Consider any triplet (^4, B, C), where A is a, p x p (complex- valued) constant matrix, 
B is apxl (complex- valued) constant matrix, and C is a 1 x p (complex-valued) constant 
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matrix. For short, we will refer to such a triplet as a triplet of size p. From such a 
triplet, construct the auxiliary pxp matrices Q and N by solving the respective Lyapunov 
equations 

QA + A^Q = C^C, 

(2.1) 

AN + NA^ = BB^. 

Note that if (Q, N) satisfies the system in (2.1), so does {Q\ N^y^ hence, there is no loss of 
generality in assuming that the solution matrices Q and N to (2.1) are selfadjoint. Then, 
form the pxp matrix F{x^ t) and the scalar quantity u{x, t) as 

F(x, t) := e2^'^-^^(^')'* + Q e'^^^-^'^'^N, (2.2) 

u{x, t) := -2B^F{x, t)-^C^. (2.3) 
Let us also define the pxp matrix G{x, t) and the scalar quantity v{x, t) as 

Gix, t) := + Are2^'^-4^(^')'*Q, (2.4) 

v{x, t) := -2CG{x, t)-^B. (2.5) 

Theorem 2.1 Given a triplet (A, 5, C) of size p, let (Q, N) be a selfadjoint solution to the 
system in (2.1). Then u{x,t) defined as in (2.3) satisfies the NLS equation (1.3) at any 
point on the xt-plane where F{x,t) is invertible. Similarly, v{x,t) defined in (2.5) satisfies 

(1.3) at any point on the xt-plane where G{x,t) is invertible. 

PROOF: Let us only give the proof for v because the proof for u is similar. In fact, v{x, t) 
turns into u{—x, t) when replacing {A, B, C) with {A\ C\ B^)., which can also be used for 
the proof related to w(a;, t). Let us drop the arguments of the functions involved and write 

(2.4) and its adjoint as 

G = e"^ + Ne^^ Q, G^ = e"^^ + Qe^iV, (2.6) 
where we have defined 

/3 := 2Ax + AiAH. (2.7) 
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Note that v in (2.5) is well defined as long as G -"^ exists. Taking appropriate partial 
derivatives, from (2.5) and (2.6) we obtain 

ivt + v^ox + 'ivvh = -2CG-^PG-^B, (2.8) 

P := -iGt - G^^ + 2G^G-^G^ + ^BB\G^)-^C^C. 

From (2.6) we get 

Gt = -AiA^e-^ - AiN{A^fel^'Q, G^, = -2Ae-^ + 2NA^ef^^Q, (2.9) 

= AA'^e-'^ + m(A^fe^^Q, (2.10) 

e^'QG-i = (G^)-iQe^, e^Ar(G^)-i = G'^Nef^'. (2.11) 

With the help of (2.1), (2.6), and (2.9)-(2.11), one can verify that P = 0, and hence the 
right hand side of (2.8) is zero. | 

There are several questions that can be raised. For example, are the Lyapunov equa- 
tions given in (2.1) solvable; if they are solvable, are they uniquely solvable? Are the 
matrices F and G defined in (2.2) and (2.4), respectively, invertible? The answers to these 
questions are affirmative under appropriate restrictions on the triplet {A, B, C), as we will 
see. 

Let us note that u and v defined in (2.3) and (2.5), respectively, are analytic functions 
of X and t at any point on the a;t-plane as long as the matrices F and G, respectively, 
are invertible at that point. This is because the entries of those matrices and hence also 
their determinants can be written as sums of products of sine, cosine, exponential, and 
polynomial functions of linear combinations of x and t. 

Consider the scalar function defined as 



n{x) := Ce-^^S. 
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(2.12) 



The right hand side is called a matrix realization of fl in terms of the triplet {A,B,C). 
Without changing fl,{x), it is possible to increase the value of p in the size of the triplet 
{A,B,C) by padding the matrices A, B, C with zeros or by modifying A, B, C in some 
other fashion (cf. [8], Subsection 2.4). Conversely, it might also be possible to reduce the 
value of p in the triplet (A, 5, C) so that the quantity Q{x) will remain unchanged. The 
matrix realization in (2.12) is said to be minimal if the value of p in the triplet {A, B, C) 
is the smallest and yet Q,{x) remains unchanged by the choice of p. The triplet (A, 5, C) 
is minimal if and only if [5] the intersections of the kernels of CA^ and of the kernels of 
B^{A^y for j = 0, 1, 2, . . . are trivial; i.e. 



e : CA^$, = for j > O} = {0} = {ry G : B\A^yr] = for j > O} . (2.13) 



It is also known [5] that the triplet yielding a minimal realization in (2.12) is unique up to 
a similarity transformation {A^ B, C) i-^ {EAE~^^ EB, CE~^) for some unique matrix E. 

The results in the next theorem are known [3] , but they are collected here in a sum- 
marized form and a brief proof is included for the benefit of the reader. 

Theorem 2.2 Assume that the triplet (A, S, C) of size p corresponds to a minimal real- 
ization in (2.12) and that the eigenvalues of A all have positive real parts. Then: 

(i) The Lyapunov equations in (2.1) are uniquely solvable. 

(ii) The solutions Q and N are p x p selfadjoint matrices. 

(iii) Q and N can be expressed in terms of the triplet {A, B, C) as 



(iv) Q and N are invertible matrices. 

(v) Any square submatrix ofQ containing the (l,l)-entry or {p,p)-entry ofQ is invertible. 
Similarly, any square submatrix of N containing the (l,l)-entry or {p,p)-entry of N 
is invertible. 




(2.14) 
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(vi) The quantities F and G defined in (2.2) and (2.4), respectively, are p x p matrices 
invertible at every point on the xt-plane. 

PROOF: By introducing the parameter a, let us write the first equation in (2.1) as 

-Q{al - v4) + {al + A^)Q = C^C, 

or equivalently as 

-{al + A^)-^Q + Q{al - A)'^ = {al + A^)-^C^C{a - A)-\ (2.15) 

where / is the px p identity matrix. Since A and {—A^) have eigenvalues on the right and 
left complex half planes, respectively, we can integrate (2.15) along a simple and positively 
oriented contour 7 lying on the right half complex plane and enclosing all eigenvalues of 
A. Thus, we obtain Q uniquely as 

Q=-^ / da(al + A^)-^C^C(aI - A)-\ 

Similarly, the solution to the second equation in (2.1) is unique and is obtained as 

N=^ [ da{aI-A)-^BB\aI + A^)-\ 

Thus, (i) is proved. From (2.1) it is seen that (Q^iV^) is a solution to (2.1) whenever 
{Q,N) is a solution and hence from the uniqueness of the solution we obtain (ii). From 
(2.14) we see that 

QA + A^Q = - ds-f e-^'^CtCe-^^ = -e-^'^CtCe-^^T = CtC, 
Jo ds I J 

where we have used the fact that all eigenvalues of A have positive real parts. A similar 
argument for N completes the proof of (iii). From their selfadjointness and positivity as 
seen from (2.14), it follows that all eigenvalues of Q and N are nonnegative. Moreover, 
(2.13) implies that zero is not an eigenvalue of Q and N. Hence Q and N are invertible. 



proving (iv) . The positivity of all eigenvalues also implies ( v) . The invertibility of F follows 
from using Theorem 4.2 of [3] in (2.2), and the proof of invertibility for G is similar. | 

The results in the next theorem are useful in extracting the scattering data for (1.2) 
from the corresponding potential u{x,t), which is also a solution to (1.1). For the benefit 
of the reader we state such results in a summarized and unified form. The proofs of these 
results are available in Theorems 3.1 and 3.3 of [3], and hence they will not be given here. 

Theorem 2.3 Assume that the triplet {A,B,C) of size p corresponds to a minimal real- 
ization in (2.12) and that the eigenvalues of A all have positive real parts. Further, assume 
that A has m distinct eigenvalues Q!i, . . . , am and the multiplicity of aj is Uj. Then: 

(i) There exists a unique triplet {A,B,C), where A is in a Jordan canonical form with 
each Jordan block containing a distinct eigenvalue and having —1 in the superdiagonal 
entries, and the entries of B consist of zeros and ones. More specifically, we have 





r^i 


.. 


- 







A2 .. 







. 


.. 


Am - 


-aj 


-1 


. 


.. - 





aj 


-1 . 


.. 








aj . 


.. 


. 





. 





B = 



B,: = 



B2 
Bm 

0- 




1 



C=[Ci C2 



-1) 



Cji Cjo. 



where Aj has size rij x Uj, Bj has size Uj x 1, Cj has size 1 x Uj, and the constants 
Cj(^nj-i) (^1^^ nonzero. 

(ii) The triplet {A,B,C) can be constructed from {A,B,C) via 



A = MAM-^, B = MSB, C = CMS, 



(2.16) 



where M is a matrix whose columns are formed from the generalized eigenvectors of 
{—A) and S is an upper triangular Toeplitz matrix commuting with A, is uniquely 
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determined by M and B, and has the form 







. 


" 









-52 . 







s = 





. 









. 


. 







'JUj 




L 








^jl 1 
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JUj -• 



for some constants djg. 

(iii) The triplets {A,B,C) and {A,B,C) yield the same solution u{x,t) to (1.1). 

(iv) The complex constants [ioij) correspond to the hound-state poles on the upper half 
complex plane of the transmission coefficients Ti and T^ appearing in (4-4)-(4-V 
Section 4- 

(v) For each j, the complex constants cjg for s = 0, 1, . . . , {nj — 1) appearing in C corre- 
spond to the bound-state norming constants associated with the bound-state pole {iaj) 
of the transmission coefficients. 

3. THE SCALAR NLS EQUATION 

Recall that we refer to (1.3) as the scalar NLS equation when is a scalar quantity and 
as the matrix NLS equation when m is a matrix quantity. In this section we exploit a certain 
symmetry in (2.2) and show that without changing the value of the scalar u{x,t) in (2.3) 
it is possible to transform the triplet (A, -B, C) in such a way that some or all eigenvalues 
of A can be reflected from the right half complex plane to the left half complex plane. The 
same result holds for the scalar v{x,t) given in (2.5); namely, it remains unchanged when 
the triplet (A, S, C) is transformed so that some or all eigenvalues of A are reflected with 
respect to the imaginary axis on the complex plane. Later, we will see that these results 
remain valid also for the matrix NLS equation when u and v are matrix valued. 

For repeated eigenvalues of A, the aforementioned transformation must be applied to 
all the multiplicities in such a way that after the transformation we should not have any 
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eigenvalue pairs symmetrically located with respect to the imaginary axis of the complex 
plane. As mentioned in Section 1, eigenvalue pairs symmetrically located with respect to 
the imaginary axis cannot yield soliton solutions to (1.1). For such pairs, it is already 
known [9] that (2.1) is not uniquely solvable. 

Let us write (2.2) as 

F{x, t) = Q \e-^^--^^^"t + Q-ie'^'^-^^(^')'*Ar-i] N. (3.1) 

Comparing (2.2) and (3.1), we next prove that u{x,t) appearing in (2.3) remains invariant 
under the transformation 

{A, B, C, Q, N) ^ {-A\ -N-^B, -CQ-\ -Q-\ -N'^), 

where all the eigenvalues of A are reflected with respect to the imaginary axis on the 
complex plane as a result of ^ i— > {—A'^). 

Theorem 3.1 Assume that the triplet {A,B,C) corresponds to a minimal realization in 
(2.12) and that all eigenvalues of A have positive real parts. Consider the transformation 



{A, B, C, Q, iV, F, G, u, v) ^ {A, B, C, Q, AT, F, G, u, v), (3.2) 

where {Q,N) corresponds to the unique solution to the Lyapunov system in (2.1), the 
quantities in {F,G,u,v) are as in (2.2)-(2.5), 

A = -A^, B = -N-^B, C = -CQ-\ Q = -Q-^, N = -N'^ (3.3) 



and (F, G, u, v) is as in (2.2)-(2.5) but by using {A, B, C, Q, N) instead of {A, B, C, Q, N) 
on the right hand sides. We then have the following: 

(a) The quantities F and G are transformed as 

F = Q-^FN-\ G = N-^GQ-\ (3.4) 
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(b) Q and N satisfy the respective Lyapunov equations 

QA + A^Q = C^C, 

AN + NA^ = BB^. 



(3.5) 



(c) The matrices Q and N are selfadjoint and invertible. 

(d) The matrices F and G are invertible at every point on the xt-plane. 

(e) u{x, t) — u{x, t) and v{x, t) — v{x, t). 

PROOF: Using (3.3) in (2.2) and (2.4) we get (a). Using (2.1) and (3.3), it can directly be 
verified that (3.5) is satisfied, proving (b). Using (ii) and (iv) of Tlieorem 2.2 in (3.3), it 
follows that (c) holds. The invertibility in (d) follows from (3.4) by using the invertibility 
of Q and N stated in (c) and the invertibility of F and G stated in (vi) of Theorem 2.2. 
Finally, (e) can be proved by using (2.3) and (2.5) with the help of (3.3), the selfadjointness 
of Q and N, and (3.4). | 

Next, we show that even if we reflect some of eigenvalues of A from the right to 
the left half complex plane. Theorem 3.1 remains valid by choosing the transformation in 
(3.2) appropriately. For this purpose, let us again start with a triplet {A, B, C) of size p 
and corresponding to a minimal realization in (2.12), where the eigenvalues of A all have 
positive real parts. Without loss of any generality, let us partition ^, S, C as 







" 






A = 






, B = 









A2 




B2_ 



c 



:ci C2], 



(3.6) 



so that the q x q block diagonal matrix Ai contains the eigenvalues that will remain 
unchanged and A2 contains the eigenvalues that will be reflected with respect to the 
imaginary axis on the complex plane, the submatrices Bi and Ci have sizes q x 1 and 
1 X g, respectively, and hence A2, B2, C2 have sizes ip — q)x{p — q), {p — q)xl, 1 x (p — q), 
respectively, for some integer q not exceeding p. Let us write the corresponding respective 
solutions to (2.1) as 



'Qi 


Q2 




"iVi 


N2 








_Q3 








iV4_ 



(3.7) 
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where Qi and Ni have sizes q x q, Q4 and A^4 have sizes {p — q) x (p — q), etc. Note that 
because of the selfadjointness of Q and N stated in Theorem 2.2, we have 

= QI^Q3, Qi = Q4, Nl^Ni, N^^Ns, iVj = iV4. (3.8) 

Furthermore, from Theorem 2.2 (v) it follows that Qi, A^i, and A^4 are all invertible. 

Let us clarify our notational choice in (3.6) and emphasize that the partitioning in 
(3.6) is not the same partitioning used in (i) of Theorem 2.3. 

Theorem 3.2 Assume that the triplet {A,B,C) partitioned as in (3.6) corresponds to 
a minimal realization in (2.12) and that all eigenvalues of A have positive real parts. 
Consider the transformation (3.2) with {A,B,C) having similar block representations as 
in (3.6), {Q,N) as in (3.7) corresponding to the unique solution to the Lyapunov system 
in (2.2), 



ii=^i, A2 = -AI Bi = B^-N2N^^B2, B2 = -N^^B2 



(3.9) 



Ci= Ci— C2Q4 ^Qsi C2 = — C2(54 ^' 



(3.10) 



and {Q, N) given as 

Qi = Qi-Q2Q4^Q3, Q2 = -Q2Qr\ = -Ql'Qs, Q4 = -^4 \ (3-ii) 

iVi = iVi - ATsAr-iATg, N2 = -N2N^\ Ns = -N^^N3, N4 = -N^\ (3.12) 



and {F,G,u,v) as in (2.2)-(2.5) hut by using {A,B,C,Q,N) instead of {A, B,C,Q, N) on 
the right hand sides. We then have the following: 



(a) The quantities F and G are transformed according to 



-Q2Q4 

-Q4"' 



/ -N2N^^ 



-N 



-1 



G 



(3.13) 



(3.14) 
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(b) Q and N satisfy the respective Lyapunov equations in (3.5). 

(c) The matrices Q and N are selfadjoint and invertible. 

(d) The matrices F and G are invertible at every point on the xt-plane. 

(e) u{x.,t) = u{x.,t) and v{x,t) = v{x,t). 

PROOF: Let us use / to denote the identity matrix not necessarily having the same 
dimension in every appearance in the proof, but that dimension will be apparent to the 
reader. Note that (3.13) and (3.14) can be verified by using (3.8)-(3.12) in (2.2) and (2.4), 
which proves (a). The proof of (b) is by direct substitution in (3.5) and by using (2.1) and 
(3.6)-(3.12) and by noting that 

(3.15) 



B = 



I -N2NI 




C = C 





'I 


-N2' 


-1 


B = 






B, 













/ 





= C 










-Qs 


-Q4 



-1 



(3.16) 



/ 

where the invertibility of Q4 and A^4 is also used, which follows from Theorem 2.2 (v). The 
selfadjointness of Q and N follows from (3.7), (3.8), (3.11), and (3.12). The invertibility 
of Q and N can be seen from (3.11) and (3.12) by writing 



Q = 



Qi Q2 
/ 



/ 

Q4^Q3 —Qa^ 



Qi ~Qi Q2 
/ 



-1 


" / " 




^Qs —Qa 



-1 



, (3.17) 



and a similar expression for N and the fact that Qi, (^4, -/Vi, N4 are all invertible. Thus, 
(c) is proved. The invertibility in (d) follows from (3.13) and (3.14) by using the invert- 
ibility of Q4 and A^4 stated in (v) of Theorem 2.2 and the invertibility of F and G stated 
in (vi) of Theorem 2.2. Finally, we evaluate F{x, t) by using {A, B, C, Q, N) instead of 
{A, B,C,Q, N) on the right hand side of (2.2) and evaluate u{x,t) by using {B,F,C) 
instead of {B,F,C) on the right hand side of (2.3). With the help of (3.7)-(3.17), it 
is straightforward to show that u{x,t) simplifies to u{x,t); similarly, v{x,t) simplifies to 
v{x,t), completing the proof of (e). | 
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Finally in this section, we show that u and v defined in (2.3) and (2.5), respectively, 
can be transformed into each other by transforming the triplet {A, B,C) in a particular 
way. 

Theorem 3.3 Assume that the triplet {A, B, C) corresponds to a minimal realization in 
(2.12) and that all eigenvalues of A have positive real parts. Consider the transformation 
(3.2) where {Q,N) corresponds to the unique solution to the Lyapunov system in (2.1) and 

{F,G,u,v) is as in (2.2)-(2.5), 

A = A, B = Q-^C\ C = B^N-\ Q = N'^ N = Q-\ (3.18) 



and (F, G, u, v) is as in (2.2)-(2.5) hut by using (A, B, C, Q, N) instead of (A, 5, C, Q, N) 
on the right hand sides. We then have the following: 

(a) Q and N satisfy the respective Lyapunov equations given in (3.5). 

(b) The matrices Q and N are selfadjoint and invertible. 

(c) u{x, t) — v{x, t) and v{x, t) — u{x, t). 

PROOF: Using (2.1) and (3.18), it can directly be verified that (3.5) is satisfied. The 
selfadjointness and invertibility of Q and follow from (3.18) because Q and have 
those properties, as stated in Theorem 2.2 (ii) and (iv); hence, (b) holds. Finally, (c) is 
proved by direct substitution in (2.3) and (2.5) with the help of the selfadjointness of Q 
and N. | 

4. THE SCATTERING COEFFICIENTS AND JOST SOLUTIONS 

In this section, we evaluate the Jost solutions and the corresponding scattering co- 
efficients for the Zakharov-Shabat system (1.2) associated with the NLS equation (1.3) 
when the potential is given by (2.3) or (2.5). In analyzing (1.2), we will use the nota- 



tion (p 



, where the subscripts 1 and 2 denote the first and second components. If 



(f{X, X, t) and (p{X, x, t) are any two solutions to (1.2), their Wronskian [ip; 0\ is independent 
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of X, where we have defined 



We stress that an overbar does not indicate complex conjugation. 



(4.1) 



Recall [13,15] that the Jost solutions i/^, i/^, 0, and (j) are defined as those vector solutions 
to (1.2) with asymptotics 



V'(A,x,t) 



0(A,a;,t) 







-i\x 







[1 + 0(1)], ^|;{X,x,t) 



[l + o(l)], ^{X,x,t) 



—iXx 





„iXx 



[1 + 0(1) 



X +00, (4.2) 



[l + o(l)], a;^-oo. (4.3) 



When tt(-, t) is integrable for each fixed t, the four Jost solutions exist and their asymptotics 
at the other end of the real axis yield the scattering coefficients R, L, Ti, via 



i){\x,t) = 
(t){\x,t) 



(T;)-ie-^^^ 
-Lt(Trt)-ie^^* 

rp-l^-iXx 



[1 + o(l)], X -oo, 

[l + o(l)], X ^ -oo, 
1 + o(l)], X — >■ +00, 



-i?+(T/)-ie-*^^ 



[l + o(l)], x^+oo. 



(4.4) 

(4.5) 
(4.6) 
(4.7) 



The scattering coefficients can equivalently be obtained with the help of the Wronskian 
defined in (4.1); in fact, we have 



LTfi = (i?r-i)t = [</,;V^], Tfi = [(^;V'], Tr^ = [V5;0]. 



(4.8) 



For the Zakharov-Shabat system (1.2) we have T\ = % if the scalar potential u{x,t) 

vanishes as a; — >^ ±oo. However, we will retain the separate notations for Ti and for a 
subsequent generalization to the matrix case, where T/ and may differ. 
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The Jost solutions ip and ip and the potential u{x,t) are recovered [13,15] as 



ip{X,x,t) 



ip{X,x,t) 



i) 

-iXx 



J X 

+ / dyK{x,y,t)e-'^y, 

J X 



u{x, t) = -2[l ] K{x, x,t) = 2 K{x, x, t)^ 

from the solutions to the Marchenko integral equations 



Qi{x + y,t) 

'ni{x + y,ty 





K{x,y,t) + 
K{x,y,t) - 



+ / dy K{x,z,t)Vt\{z + y,t) ^ 
dy K{x, z, t) Qi{z + y, t)^ = 



(4.9) 

(4.10) 
(4.11) 



y>x, (4.12) 
y>x. (4.13) 



The Jost solutions ip and ip corresponding to u{x,t) given in (2.3) with [A^B^C) 
(^i,5i,Ci) can be evaluated by solving (4.12) and (4.13) with 



Since 



l]i(y,t) = C,e-^'^-4^^?*5i. 



ni{x + y,t) = Cie^^'^e-^'^-^^^' 



(4.14) 



ipi^X, x, t) 



ip{X,x,t) 



the Marchenko equations in (4.12) and (4.13) have integral kernels separable in z and y 
and hence they can be solved explicitly by algebraic methods, yielding 

iBlFi{x,t)-\XI + iAl)-^e'^''Cl 

e»^^-zCi[Fi(a;,t)t]-iArie-2A^^+4i(A^)'t(A/ + af)-ie^^^C;. 

>-^^- +iB{e-^^l''+^'^^l^''Qi[F^{x,ty]-^XI ~ iAiy^e-'^^'Bi 
iCi[Fi{x, t)^-^{XI - ai)-ie-*^^Si 

Similarly, the Jost solutions (p and and the potential v{x,t), in the form given in 
(2.5), are recovered as 



(4.15) 



(4.16) 



(p{X,x,t) = 



" g — i Ax " 











J — oo 
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-i\y 



(4.17) 



4>{X,x,t) 





-,i\x 



+ j dyM{x,y,t)^ 



My 



v{x, t) = 2[l ] M{x, X, t) = -2 M{x, x, t)^ 
from the solutions to the Marchenko integral equations 



M{x,y,t) + 



M{x,y,t)- 



n,.{x + y,t) 




n,{x + y,ty 



+ dyM{x,z,t)nr{z + y,t) = 



(4.18) 
(4.19) 



- / dyM{x,z,t)n,{z + y,t)^ 

J — oo 



y<x, (4.20) 



y <x. (4.21) 



Corresponding to the potential v{x, t) given in (2.5) with {A, B, C) = (Ar, -Br, C-c) we have 



(4.22) 



(t>{X,x,t) = 



(4.23) 



(4.24) 



which yields an integral kernel separable in z and y in (4.20) and (4.21) and hence allows 
us to solve those integral equations explicitly by algebraic methods, yielding 

iB^,[Gr{x,t)^]-\XI + iAl)-^e-'^^Cl 

iC,G,{x,t)-\XI -iA^y^e'^^'B, 
_e^^^ + iSt[G'r(a;,t)t]-igre2^--*+^^^'*(A/ - iAry^e'^'^Br 

Proposition 4.1 Assume that the triplet {A,B,C) corresponds to a minimal realization 
in (2.12) and that all eigenvalues of A have positive real parts. Let {Q,N) correspond to 
the unique solution to the Lyapunov system (2.1) and let {F,G) be as in (2.2) and (2.4). 
We then have: 

(a) F{x,t)~^ and G{x,t)~^ as x ^ ±oo. 

(b) e-2^'^+4^(^')'*g[F(a;,t)t]-i ^ iV-i as 

(c) [F(a;,t)t]-iAre-2^^^+4»(^^)'*^g-i as 

(d) [G(a;,t)t]-iQe2^^+4^^'* ^ N'^ as x ^ +oo 
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X —>■ — oo. 



X — > — oo. 



(e) G'(a;,t)-iiVe2^^^-4»(^^)'* ^ Q-i as x ^ +00. 

PROOF: The proofs involving G are obtained by using (2.6) the same way as in the proofs 
for F, and hence we will only give the proofs for F. Using (2.2) and (2.7) we see that 

F = V Qe-^N, F^ = + Ne'^^ Q, 

and hence 

F = ef^\l + e-f^^Qe-f^N], Q'^ FN'^ef^ = I + Q-^e^^^ N'^e^^ , (4.25) 

F^Q-^ef^^N-^ = [I + ef^Q-^ef^^N-^], e^^ N'^F^ = [I + e^^ N-^e^Q-^]Q. (4.26) 
Forming the inverses from (4.25) and (4.26) we obtain 

F-^ = [I + e-l^'Qe-^N]-^e-^\ F'^ = N'^e^il + Q'^e^' N-^e^]-^Q-\ (4.27) 

e-^'Q(Ft)-^ = N-\I + e^Q-^e^'N-Y\ (4-28) 

(Ft)-iiVe-^^ = Q-^[I + e^^ N-^e^Q-^]~^. (4.29) 

By letting x — > ±00 in (4.27) we prove (a). By letting x —00 in (4.28) and (4.29) we 
prove (b) and (c), respectively. The remaining proofs involving G are obtained in a similar 
manner. | 

The following result is known [5], but it will be useful later and hence we provide a 
brief proof. 

Proposition 4.2 For any triplet of matrices (A, S, C) with sizes p x p, p x n, and n x p, 
respectively, we have 

[In + C{XIp - A)-^B]-^ =In- C{XIp -A + BC)-^B. (4.30) 

PROOF: Using the identity 

EC = (XIp -A + BC) - (XIp - A), 
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one can verify that the product of one side of (4.30) and its inverse yields /„. | 

For convenience we also give a short proof of the following already known result on 
Schur complements [9]. 

Proposition 4.3 For any matrices U and V with sizes n x p and p x n, respectively, we 
have the determinant identity 

det(/„ + UV) = det(/p + VU). (4.31) 

PROOF: Let us use the decompositions 
'ip + VU 

_Onp lu + UV 

where O^fc denotes the zero matrix of size j x k. Note that the right hand sides in (4.32) 
and (4.33) have the same determinant, and by equating the determinants of the left hand 
sides we get (4.31). | 

Theorem 4.4 Assume that the triplet (A, B, C) corresponds to a minimal realization in 
(2.12) and that all eigenvalues of A have positive real parts. Let (Q, A^) correspond to the 
unique solution to the Lyapunov system in (2.1), {F,u) be as in (2.2) and (2.3), and the 
scattering coefficients he defined as in (4.4)-(4-V- then have: 

(a) u{x,t) —>^0 as X ^ ±oo. 

(b) For A e R, the transmission coefficients T\ and T^. appearing in (4-4) osnd (4-6), 
respectively, and their inverses are given by 

Tf ^ = 1 - iCQ-^{XI + iA^)-^C^, (4.34) 

Ti = l + iC{XI -iA)-^Q-^C\ (4.35) 
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Ip -V 



Ip V 
-u /„, 



Ip 

u 



pn 



(4.32) 



Ip ^pn 

U h. 



Tp V 

-u h 



Ip -V 

Onp In 



(4.33) 




l-iB^(XI + iA^)-^N-^B, 



(4.36) 



Tr = l + iB^N-\XI -iA)-'^B, 



(4.37) 



and hence they are functions of A alone and do not depend on t. 

(c) The reflection coefficients L{\,t) and R{X,t) appearing in (4-4) '^n'd (4-6) are both 
identically zero. 

(d) The transmission coefficients can he written as the ratio of two determinants as 



PROOF: Let us compare (4.4) and (4.5) with (4.15) and (4.16) by ignoring the subscript 1 
in (4.15) and (4.16). Using Proposition 4.1 (a) in the first component of (4.15), wc sec that 
LT^^ = 0. Then, with the help of (4.8) we also conclude RT~^ = 0. Using Proposition 4.1 
(c) in the second component of (4.15) we obtain (4.34), and similarly by using Proposition 
4.1 (b) in the first component in (4.16) we obtain 



Applying Proposition 4.2 on (4.34) and (4.36) and simplifying the resulting expressions 
with the help of (2.1), we obtain (4.35) and (4.37). We then get (c) from LTf ^ = and 
RT~^ = 0. Applying Proposition 4.3 onto (4.35) and (4.37) and using (2.1) we obtain 
(4.38). I 

5. THE MATRIX NLS EQUATION 

In this section we show that all the results presented in Sections 2-4 for the scalar 
NLS equation remain valid for the matrix NLS equation as well. In other words, when the 
scalar function u{x,t) in (1.3) is replaced with an m x n matrix-valued function of x and 
t, our results in the previous sections remain valid. 



det(A/ + zylt) 
det(A/ - lA) ' 



rr(A) 



det(A/ + at) 
det(A/ - lA) ' 



(4.38) 



{T-y = 1 + iB^N-\XI - iA)-^B. 



(4.39) 
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In Section 2 we have started with the triplet {A,B,C), where the sizes of A, B, C 
were p x p, p x 1, and I x p, respectively. In Sections 2-4 we have carefully stated all our 
results so that they all will remain valid also for the m x n matrix NLS equation in (1.3). 
In order for the reader to see the validity of all the results in Sections 2-4 in the matrix 
case, in this section we will use subscripts to indicate matrix dimensions by writing A as 
App, B as Bpm, C as Cnp, etc. At appropriate places by interpreting 1 as the identity 
matrix 1^ or and interpreting as one of the four zero matrices 0mm j Onn; 0^^, 0„^, 
we will show that all the results in Sections 2-4 remain valid also in the matrix case. 

To produce exact solutions to the m x n matrix NLS equation (1.3), we start with a 
triplet {App, Bpm, Cnp)- The Lyapunov equations in (2.1) remain unchanged with Qpp and 
Npp as the solutions. Theorem 2.1 and its proof remain unchanged, with the matrices Fpp 
and Umn still defined as in (2.2) and (2.3), respectively. Note that v{x,t) defined in (2.5) 
now becomes an n x m matrix, and Theorem 2.1 remains valid with Vnm as in (2.5) and 
Gpp as given in (2.4) for the n x m matrix NLS equation 

dv d'^v + 

Theorem 2.2 and its proof hold verbatim with Qpp, Npp, Fpp, Gpp in the matrix case. 
Alternatively, in defining v one can simply use a triplet {A, B, C) with sizes p x p, p x n, 
m X p, respectively, so that v has the same size as u, in which case Theorem 2.1 and 
its proof hold verbatim for the m x n matrix NLS equation. In fact, this alternative is 
preferred, because in that case the generalizations of all the results in Sections 2-4 to the 
matrix case are simpler to state. Theorem 2.3 has a counterpart in the matrix case, but 
we will not consider that generalization in this paper. 

In terms of Umn, Vnm, Umn, Sixid Vnm, cih the results in Section 3 until Theorem 3.3 also 
hold verbatim in the matrix case. Theorem 3.3 and its proof also hold when stated with 
iJ'mnj Vnm J '^nmj and Vmn- Alternatively, in evaluating v and v, by using a triplet (^4, B, C) 
with sizes p x p, p x n, m x p, respectively, all the results in Section 3, including Theorem 
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3.3 and its proof, hold verbatim in the matrix case. 



All the results in Section 4 also hold with obvious interpretations of the stated quanti- 
ties. The matrix Zakharov-Shabat system corresponding to (1.2) now becomes (1.4). The 
Jost solutions ^A(m+n)n, <P{m+n)m, i^{m.+n)7n, (p{m+n)n are uow matrices. In analogy with 
(4.2) and (4.3), their asymptotics in the matrix case are given by 
0. 



(f){X,x,t) 



pi\x T 
„—i\x T 



Or 



[l + o(l)], i;iX,x,t) 



[l + o(l)], ^{X,x,t) 



„—iXx J 
c J-m 



0. 



Omn 

„i\x J 
e -In 



[l + o(l)], X^+OO, 



[l + o(l)], x^-oo. 



The definition of the Wronskian in (4.1) holds in the matrix case, where it is now a square 
matrix of size either m x m or n x n, depending on the sizes of matrix solutions. The 
asymptotic relations (4.4)- (4. 7) defining the scattering coefficients hold as stated, where 
now the scattering coefficients Rnmj (^i)nn) and {Tj.)mm are matrices. The Fourier 

relations (4.9) and (4.10) are also valid in the matrix case and are given by 

OrrjTj, 



ip{X,x,t) 



ilj{X,x,t) 



„i\x f 



-iXx ; 



+ / dyK{x,y,t)e'^y, 



Or 



+ / dyK{x,y,t)e 



-iXy 



where K{x, y, t) is now an {m + n) x n matrix and K{x, y, t) is an {m + n) x m matrix. 
The Marchenko equations in (4.12) and (4.13) remain unchanged with the understanding 
that the quantity fii(y, t) appearing in the kernel is now an n x m matrix and we have 



K{x,y,t) + 
K{x,y,t) - 



0mm 

Qi{x + y,t) 
'ni{x + y,t)^ 

Onn 

and the potential Umn is recovered as in (4.11) via 



dy K{x, z, t) Qi{z + y,t) = 
dyK{x, z,t) ^li{z + y^t)"^ 



On 
Or 



y> X, 



y> X, 



u{x, t) = -2 [ Im Omn ] K{x, X,t) ^ 2 K{x, X, t)^ 
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Similarly, the Fourier representations in (4.17) and (4.18) remain true by interpreting 
them as 



(l){X,x,t) 



(/!)(A, X, t) 







nm 





J — oo 



mn 



„iXx T 
c -In 



+ / dyM{x,y,t)e 



iXy 



where now M{x, y, t) and M(a;, y, t) have sizes (m + n) x m and (m + n) x n, respectively. 

In generalizing (4.22) to the matrix case, let us simply use a triplet {Aj-^Br^Cj.) with 
sizes pxp, pxn, mxp, respectively. With the understanding that t) is now an m x n 
matrix, the Marchenko equations (4.20) and (4.21) hold true and can also be written as 



M{x,y,t) + 



M{x,y,t)- 



fl^{x + y,t) 

Omm. 



+ / dyM{x,z,t)nr{z + y,t) 



dyM{x,z,t) Qr{z + y,t)^ 



On, 

0. 



y <x, 



y <x, 



Qr{x + y,ty 

and the potential Vmn is recovered as in (4.19) via 

V{x,t) = 2[Im Omn]M{x,X,t) = -2M{x,X,t)^ 



When the triplet {A\, B\, C\) with matrix sizes p x p, p x m, n x p, respectively, is used 
in the Marchenko kernel (4.14), the Jost solutions '4'{m+n)n ^tnd '4^{rn+n)m on any half line 
X e (a, +oo) for any real number a can be expressed as in (4.15) and (4.16) via 

zS/Fi(a;, t)-i(A/ + af)-ie^^^C/ 
.e»^^/n - iCx[Fx{x, t)t]-iiVie-2^i^^+^*(A^)'*(A/ + af)-ie»^^C/ , 

■e-^^^/„, + z5;e-2A^^+4»W)'*Qi[Fi(a;,t)t]-i(A/-ai)-ie-^^^Si" 

Similarly, when the matrix triplet (Ar, Bj.,Cj-) of sizes pxp, pxn, mxp, respectively, 
is used in the Marchenko kernel (4.22), the corresponding Jost solutions 0(m+n)m and 

24 



^(A,x,t) 



ip{X,x,t) = 



(f)(rn+n)n On any half line x e (—00, a) for any real number a can be expressed as in (4.23) 
and (4.24) via 

'e-^^^^I^ -zCrG'r(a;,t)-iArre2^^=^-4*(^^)'*(A/ + at)-ie-^^^C't' 
z5t[Gr(a;,t)t]-i(A/ + at)-ie-^^^Ct 



4>{\,x,t) 



(f){X,x,t) 



iCrGr{x,t)-HXI -iA, 



-1 „iAa 



Proposition 4.1 and its proof hold as stated with the understanding that we interpret 
in (a) as 0pp. Proposition 4.2 already contains the matrix case in (4.30), and Proposition 4.3 
is already stated in such a way that it can directly be used in the matrix case. Theorem 4.4 
and its proof in the matrix case hold with the following minor diflFerences. The quantities 
u, Ti, Tr, L, R are now matrices Umn, {T\)nn, {Tr)mm, Lmn, Rnm, respectively; (4.34)-(4.37) 
and (4.39) hold as stated with the understanding that the scalar 1 becomes an identity 
matrix / of appropriate size, namely, 

Tf 1 = - iCQ-\XIp + iA^)-^C^, 

Ti = In + iC{XIp - iA)-^Q-'^C^, 

= Im- iB^XIp + iA^)-^N-^B, 

T, = Im + iB^N-\XIp - iA)-^B. 

In the matrix case, on the left hand sides of the two equations in (4.38) we need to replace 
Ti and Tj. by their determinants; namely, (4.38) in the matrix case becomes 

det r,(A) = + det r,(A) = + '^*> 



det(A/p - iA) 



det {XI p - iA) 



6. AN EXAMPLE 

We conclude our paper by providing an application of Theorem 3.2 to a specific case. 
In Example 7.2 of [3] we evaluated the exact solution to the NLS equation corresponding 
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to the triplet {A, B, C) with 





"2 " 




"l" 


A = 




, B = 






-1 




1 



c 



(6.1) 



We incorrectly conjectured in that example that we had a nonsoliton solution because one 
of the eigenvalues, of A was not positive. Using the result of Theorem 3.2 of the present 
paper, we are now able to confirm that that solution is indeed a two-soliton solution. For 
this, we proceed as follows. 

Using the triplet {A^B^C) of (6.1), we solve (2.1) in a straightforward manner and 

get 

Q = 



"1/4 -1 


, N = 


"1/4 1 


-1 -1/2 




1 -1/2 



We then construct (A, S, C, Q, TV, F) via (3.9)-(3.13) and obtain 

C'=[3 -2], 





"2 


0" 




"3" 


i = 






, B = 









1 




2 



(6.2) 



Q = 



"9/4 -2" 




'9/4 2" 




, N = 


-2 2 




2 2 



F = 



Ax-im I ^ p 
^16 



-4x-16it 



9 



4e 



p _j_ ^p—2x—Ait Jlx—Ait _ ^g— 4a; — 16it i ^p — ix—Ait 

2" 



-2x-Ait 



-Ax-l&it 



4e 



-2x-Ait 



As seen from (6.2), the eigenvalue (—1) of A is transformed into the eigenvalue (+1) of A. 
The potential u{x,t), or equivalently u(x,t), is then constructed via (2.3) and we get 



Ait I 



u{x, t) 



(9e-^^ + 166*^) - 32ei^**(4e-^^ + 9e 
-128 cos(12t) + 4e-6^ + 166^^ + 81e-2^ + 64e2^ 



agreeing with the potential of Example 7.2 of [3]. For this potential, the transmission 
coefficients are evaluated via (4.38) as 



Ti(A)=Tr(A) = 



(A + 2z)(A + z) 
(A-2z)(A-z)' 



26 



because the real parts of all eigenvalues of A in the associated triplet {A, B, C) are positive. 

As for the norming constants associated with the bound states at A = 2i and A = i, we 

need to transform the triplet {A,B,C) of (6.2) into another triplet {A,B,C), different 

1' 



from (6.1), so that we will have B 



1 



By using (2.16) we obtain 





"l 


o" 




"3 


o" 




"2 


o" 




"l" 


M = 


















, B = 









1 







2 







1 




1 



C=[9 -4]. 

Thus, the norming constant at A = 2z is 9 and the norming constant at A = z is —4. Finally, 
the Jost solutions i/ji^j ^i t) ^tnd '(/'(A, x, t) to the Zakharov-Shabat system (1.2) are obtained 
via (4.15) and (4.16), respectively, by using {A, B,C,Q, N, F) for (ii, A, C*!, Qi, M, Fi) 
there. For example, for the Jost solution ip we get 

■4ze-4^+4^*£/i(A,x,t)' 
4i g2{\ x,t) 







+ 



(A + 2i)(A + i)[-128cos(12t) +4e-6^ + 166^^ + 81e-2^ + 64e2^] ' 



where we have defined 

^i(A, X, t) := 36(A + i)e^^+^2** + 16(A - i)e'''+^'^'' - 16(A + 2i)e^^ - 9(A - 2z), 
^2(A, X, t) := 48(A + i)e^'^'^ + 48(A + 2i)e-^2it _ sAe"^^ - 81(A + i)e-2^ - 32(A + 2^)6^^. 
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